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Abstract 

We study the radiative decay properties of the charmonium-like X, Y and Z mesons generated 
dynamically from vector meson-vector meson interaction in the framework of a unitarized hidden- 
gauge formalism. In the present work we calculate the one- and two-photon decay widths of the 
hidden-charm Y(3940), Z(3930) (or X(3915)) and X(4160) mesons in the framework of the vector 
meson dominance formalism. We obtain good agreement with experiment in case of the two photon 
width of the X(3915) which we associate with the 2 + resonance that we find at 3922 MeV. 

PACS numbers: 12.39.Fe, 12.40.Vv, 12.40.Yx, 13.75.Lb, 14.40.Rt 
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I. INTRODUCTION 



With the discovery of many new unexpected charmonium-like resonances in the last 4-5 
years charmonium spectroscopy got back in the focus of interest. In the past the understand- 
ing of meson structure was to a large extend based on the constituent quark model. However, 
it turned out that the meson spectrum is much richer than one might expect from the simple 
quark model predictions. This fact was first observed in the light meson sector but later on 
was also found in the heavy quark sector e.g. charmonium spectrum. The majority of the 
'new' charmonium-like X, Y and Z mesons mainly discovered at the B-factories BELLE and 
BaBar cannot be easily accommodated in the qq model and are therefore interesting objects 
for meson structure beyond the constituent quark model. The structure assumptions which 
have been studied range from threshold effects and hybrids to tetraquarks, meson-meson 
bound objects and dynamically generated resonances from meson-meson interaction. An 
overview of the present situation is e.g. given in 

In this context, the coupled channel approach which combines chiral dynamics with a 
unitarization formalism [2| turned out to provide a useful tool to determine the mass and 
width of resonances as well as information on their decay patterns, cross sections etc. In 
the coupled channel formalism resonances are generated dynamically from meson-meson or 
meson-baryon interaction and appear as poles in the corresponding scattering amplitude. 
This approach was first successfully used for the description of pseudoscalar-pseudoscalar 
meson interaction and pseudoscalars interacting with baryons which therefore suggested a 
extension to vector meson interaction. Since the chiral Lagrangian only considers interaction 
with pseudoscalar mesons, the coupled channel approach has to be extended accordingly. 
The hidden gauge formalism 3 (| provides a consistent method to include vector meson 
interaction in the above mentioned chiral unitary approach to meson-meson interaction. 
This method was first applied to pp meson interaction, leading to dynamically generated pp 
resonances around 1.3-1.5 GeV which could be assigned to the /o(1370) and /2(1270) |7|. 
Subsequently the formalism was extended to SU(4) in order to study the interesting hidden 
charm resonances around 4 GeV. In [8] the coupled channel study led to 5 resonances where 



three of them are good candidates for the Y(3940), Z(3930) and X(4160) mesons discovered 
by BELLE and BaBar js-12]. Two further predictions for resonances have not been seen so 
far and might be objects for future experimental research. 
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A further interesting topic is the study of radiative decay properties which are also a 
crucial test for hadron structure. In the hidden gauge formalism the electromagnetic in- 
teraction is included by using vector meson dominance (VMD) that is the photon couples 
to the resonance via the p,U),<fi or J ftp vector mesons in the respective coupled channels 



13|4l5|. 



In the present paper we concentrate on the hidden-charm resonances around 4 GeV 
analyzed in [8] and study the two-photon and photon-vector meson decay properties. The 
outline is as follows: In the next section we introduce the chiral unitarity approach in 
combination with the vector meson dominance formalism we use to study radiative decays 
of dynamically generated resonances. In section II III we present our results for the 77 and 
V7 radiative decays of the charmonium-like F(3940), Z(3930) and X(4160) mesons. We 
compare the decay properties with available data and other theoretical approaches provided 
they exist. Finally we summarize our work in section HVl 

II. COUPLED CHANNEL APPROACH 

In the coupled channel model resonances are generated by meson-meson interaction under 
consideration of unitarization [2]. The scattering amplitude T is set up by means of the 
Bethe-Salpeter equation 

T=(1-VG)~ 1 V, (1) 

where the kernel V is provided by the interaction between different meson channels (see Fig. 
[I]) and G is a diagonal matrix of the respective meson loops. The meson loop integrals are 
UV divergent and are regularized by either cutoff or dimensional regularization. For further 
detail we refer to 0, S]. In the present case we consider coupled channels between vector 
mesons, where the interaction is given by the hidden gauge Lagrangian. We include four 
vector contact interaction given by 

4// = y(W^ - WT) , (2) 

which is diagrammatically represented by Fig. (Ha). is the SU(4) matrix containing the 
interacting vector mesons. Further on we consider vector meson exchange terms (see Fig. [T] 
b)) which are based on the three vector meson interaction Lagrangian 

CfP = ig((d^V v - dvVjWV) . (3) 
3 



These two mechanisms are responsible for the generation of resonances as indicated in Fig. 
[2] provided that the interaction is sufficiently strong. 
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FIG. 1: VV interaction diagrams. 

While four vector contact terms and vector-exchange terms determine the mass and 
part of the width of the resonances, the two pseudoscalar decay modes are only relevant 
for the generation of the widths of the resonances. The corresponding box diagrams with 
intermediate pseudoscalars in Fig. [1] c) are based on the interaction Lagrangian 



C v ^ = -ig{V^,d^]). 



(4) 



The transition amplitude T between the initial and final coupled channels via resonance R 

V V 
R 



V 



FIG. 2: Dynamically generated resonance from Bethe Salpeter equation 



from the Bethe Salpeter equation in Eq. ([!]) can be approximated close to a pole by 



where 
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(5) 

(6) 

(7) 
(8) 



denote the spin projection operators of the respective coupled channels. Here e(l) and e(2) 
are the polarization vectors of particle 1 and 2, respectively. The indices i,j and m run over 
the spatial coordinates, i.e. m = 1, 2, 3. 



The couplings of the resonance R to each particular VV channel can be extracted from the 
residues of the T-matrix at the poles because of the approximation in Eq. (JSJ) • The coupling 
strength to each particular channel gives an idea of the importance of its contribution to 
the resonance. In case of the F(3940) and Z(3930) the D*D* coupling is dominant, which 
implies, e.g., that the F(3940) is dominantly a D*D* bound state, in agreement with the 
D*D* molecular interpretations in Q, Q while the X(4160) is mainly D* S D* S . 

We study the radiative decays of dynamically generated resonances by coupling the pho- 
ton via intermediate vector mesons as depicted in Fig. |3j The resulting amplitude of the 
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FIG. 3: Photon coupling via VMD 



I/7 transition is provided by the vector meson dominance formalism which leads to 



t Vl = Cv/ 7 -M^e M (I/)e M (7) with C Vj 



7=2 ioT P° 

^2 ioTUJ 
— ~ for 



(9) 



2 

< 3 



for J/ ip . 



where we fix g by g = ^j- with f T = 93 MeV. In case of the charmonium, i.e. J/ip, we 
consider SU(4) breaking effects by using g = g Vc = Mj/^/(2f Vc ), where f Vc = 420/a/2 MeV 



is taken from 



18|. 



The amplitude for the one- and two-photon decays are therefore given by 



Vi,V 2 



T^HV^) ex ^^) a (-^)/ lr ■ !■) 



(10) 
(11) 



where gv^ is the coupling of the resonance R to the V1V2 channel (see Tables HI and Hill and 
Fj represents the respective isospin Clebsch Gordan coefficients of the V1V2 component for 
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a certain isospin state. The isospin states are given by 

— (\D* + D*- 



L) L) , 1 = 
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\ppJ = 


o,/ 3 


= o> 
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= o> 



(12) 



(14) 



—=(\D* + D*-) - \D*°D* )) , (13) 

-^(|pV> + |p-p + ) + |pV», 

-^(|PV>-|P-P + )), (15) 

when using the phase convention (— D*°, D* + ), (D*~ , Z)*°) and (— p + , p°, p~) for the respec- 
tive isospin doublets and triplet. After summing over the intermediate vector polarizations 
in Eqs. (fTUj) and ( ITT]) , the amplitudes for the i? — > 77 and i? — >■ V^7 decays are finally given 
by 

43 = f E ^^?A 7 x Ff x F y , , (16) 

T i? = 3 E 0^P#<W^7 x Fj x F 77 , (17) 

^ Vi,y 2 = i o ,w,(j>,j/ii> 

where Pvvt and vffi are the spin projection operators of Eqs. © and (|HJ) with the 
polarizations of V7 in the first case and of 77 in the second case. Due to the use of the 
unitarity renormalization with an extra factor 4g in case of identical particles we need 
to correct this factor when calculating observables. The unitarity normalization and the 
symmetry factors are combined in F 77 and i 7 V 7 with 

, y/2 for a pair of identical particles, e.g p°p° 
F 77 = { P y (18) 

2 for a pair of different particles, e.g p°u , 

\/2 for a pair of identical particles, e.g p°p° 
'Vj = S (19) 
1 for a pair of different particles, e.g p 00 . 

The radiative decay widths T 77 and IV 7 can be easily calculated from the transition ampli- 
tudes T by the relations 

^ = 25+116^2 ' ^2 l T 77 } | > ( 2 °) 
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where the summation over all spin states contributes the factors 
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III. RESULTS 



In order to compute the decay widths of the F(3940), Z(3930), X(4160) and the so far 
not observed '1^(3912)' we take the couplings between resonances and the VV channels 
from [8]. The couplings for isosinglet resonances are given in Tab. [J while the couplings of 
the 1 = 1 state are indicated in Tab. [Til In the case of the state with quantum numbers 
I G (J PC ) = CT (l + ~) predicted in |8j all couplings to vector mesons with hidden flavor are 
zero due to C-parity violation. Therefore radiative decays via VMD are forbidden in this 
case. 

The results for the radiative decay widths are summarized in Tab. IIHI The decay 
widths of the Y(3940) are in general smaller compared to the other resonances. A fur- 
ther common feature is that the p7 and 77 decay modes are suppressed in comparison to 
the to>7 and 07 decays except for the predicted 1^,(3912) resonance, which shows a rather 
strong coupling to the p7 decay channel. The relatively small two-photon decay width 
T(y(3940) — > 77) = 0.013 keV underestimates the corresponding result in the hadronic 
molecule interpretation of r(F(3940) — > 77) = 0.33 keV for J PC = ++ in 17| by more 
than one order of magnitude. However, the coupled channel analysis also regards the D*D* 
component as dominant. Yet, we observe that in terms of VV loops to which the two 
photons couple, there is a strong cancellation between the contribution of D*D* loops and 
pp loops, the latter ones not considered in 13]. Even when restricting to the D*D* loop, 
substantial differences from the results of jl7] remain. As a consequence of the cancellations 
found in the 77 rate, the uncertainties can be regarded to be large, we deal with about a 
factor two or three. The differences with the results of [17]] might be an indication of proba- 
bly inaccurate approximations in our scheme (we will improve this result below). One of the 
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pole positions and 7 [J ] 






3943 + i7.4, 0+[0++] 


3922 + i26, 0+[2++] 


4169 + i66, 0+[2++] 


PP 


-22 + i47 


-75 + i37 


70 + i20 


0000 


1348 + i234 


1558 + il821 


3 - i2441 


# 


-1000 - il50 


-904 - il783 


1257 + i2866 


WW 


417 + i64 


1783 + H97 


2681 + i940 


00(f) 


-215 - il07 


91 - i784 


1012 + il522 


UjJ/lp 


-1429 — *216 


-2558 - i2289 


-866 + «2752 


<pj/i; 


889 + »196 


918 + i2921 


-2617-t5151 



TABLE I: Couplings gt in units of MeV for the resonances with 1 = 0. 



pole position 


I G [J PC ] 


PP 


poo 




p4> 


3919 + i74 


1"[2++] 





-1150-z3470 


2105 + i5978 


-1067-i2514 



TABLE II: Couplings g { in units of MeV for I = 1, J = 2. 

approximations done in 8] is to neglect the three momenta of the particles with respect to 
the vector mass, \p\/Md* ~ 0. This is indeed the case for massive states but not if we deal 
with two photons in the final state since \pj\ ~ Md*. Hence, improvements on the present 
results are necessary. For this purpose we reanalyze the couplings of js| by considering the 
final momenta of the light vectors which finally couple to the photons. We therefore modify 



the corresponding expressions in js] in the following way. suited for the use of the Coulomb 
gauge for the photons. 
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pole [MeV] 


jG jPC 


meson 


r p7 [keV] 


iWkeV] 


I> 7 [keV] 


r J/l/)7 [keV] 


T 77 [keV] 


(3943, +i7.4) 


0+ (0++) 


F(3940) 


0.015 


0.989 


13.629 


0.722 


0.013 


(3922, +i26) 


0+ (2++) 


Z(3930) 


0.040 


15.155 


95.647 


13.952 


0.083 


(4169, +i66) 


0+ (2++) 


X(4160) 


0.029 


10.659 


268.854 


125.529 


0.363 


(3919, +i74) 


1- (2++) 


%(3912)' 


201.458 


114.561 


62.091 


135.479 


0.774 



TABLE III: Pole positions and radiative decay widths. 

1) For J = 0, we use: 

(*? + fc°)(* a + k°) -> (hi + k° 3 )(k° 2 + kl) + hkl + (24) 

exc 

2) For J = 2, we change: 



a 1 /-! /^X/^-i A' 1 , 1 ) — \j o? + -jrkl (25) 



with 



a = (k° 1 + k°)(k° 2 + k'l) + kl + ^, (26) 

exc 

where A; 3 is the photon momentum and M exc is the mass of the exchanged vector 
meson, Mp* or Mr>*. 

With this prescription we run the coupled channels Bethe-Salpeter equations and obtain 
new couplings of the resonances to pairs of vectors, which, via the VMD, provide the ap- 
propriate couplings of the resonances to the photons. The couplings to D*D* and D*D* 
are practically unchanged, as well as the masses of the resonances. However we obtain new 
effective couplings to the light vectors which characterize the 77 decay rates. The new re- 
sults are given in Table IIVI Due to the inclusion of the three-momenta of the photons the 
decay widths undergo changes of 30% in all cases, except for the Y(3940) state, where the 
abnormally small former width becomes by a factor of six larger. Our new result lies much 
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pole [MeV] I G J PC meson r° 7 w [keV] 



(3943, +i7A) 




- (0++) 


y(3940) 


0.085 


(3922, +i26) 


0+ 


- (2++) 


Z(3930) 


0.074 


(4169, +i66) 


0+ 


- (2++) 


X(4160) 


0.54 


(3919, +i74) 


1" 


(2 ++ ) 


%(3912)' 


1.11 



TABLE IV: New values of the r 77 w after running the coupled channel Bethe-Salpeter equations 
with the modifications of Eqs. (|25p and (|26[) . 



closer to the prediction of |17j based on a D*D* molecular structure interpretation of the 
F(3940). The two photon width obtained within the coupled channel formalism is still by 



a factor three smaller than the result of the effective Lagrangian approach of 17|] but this 
discrepancy lies within the theoretical uncertainties of this small rate. The inclusion of the 
photon three-momenta also leads to an increase of the jV decay widths quoted in Table II III 
by about a factor of two in case of the Y(3940). However this deviation can be absorbed in 
the uncertainty of our results. 

For the two-photon width of the X(4160) we obtain = 0.54 keV. In the present 
coupled channel approach the X(4160) is found to be dominantly a D*D* state. This is 



the same underlying structure as the D* S D* S bound state studied in [17j. In reference 17] 



;he D*D* molecular state was associated with the narrow Y(4140) discovered by the CDF 
19| because it was possible to explain the sizable observed J/ip(f) decay width of this state. 
Our association to the broader X(4160) is suggested by the large total theoretical width 
which was not evaluated in 17]]. Since the nature of the resonances is the same in both 



approaches it is less surprising that our result for T 77 = 0.54 keV agrees with the T 77 = 0.5 



keV evaluated in 



Experimental observations concerning the radiative decays are rare. However, the BELLE 
Collaboration searched for charmonium-like resonances in the 77 — > uJ/ip process (20! which 
resulted in an enhancement of the cross section around M = 3915 ± 3 ± 2 MeV. The peak 
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was associated to a resonance denoted by X(3915). But it is thought that it could be the 
V(3940) resonance, or even the Z(3930) which we have associated to our J p = 2 + resonance 



at 3922 MeV. In 



20[, the X(3915) has unknown spin and parity, but + or 2 + are preferred. 



In the following we compare the ratios 

' (61 ± 17 ±8) eV for J p = 0+ 
r 77 (X(3915))£(X(3915) uJ/iJ>) = { (27) 

(18 ± 5 ± 2) eV for J p = 2+ 



quoted in [20| with the results of the present approach. Let us evaluate Eq. ( 1271) for the two 
theoretical states 0+ at 3943 MeV and 2+ at 3922 MeV in Table ED By using the simple 
formula for the decay width to ojJ/ip: 



1 kg t 



2 

r ^ ~ ' (28) 

with k the momentum of the final meson, and taking the couplings from Table HJ we 

obtain 

r w j/^(o+,3943) = 1-52 MeV 

r UtW 2+,3922) = 8.66 MeV. (29) 
Together with the two photon decay widths of Table HVl we find 

r 7r B((0 + ,3943) ^cuJ/ip) = 7.6 eV 

r 7r B((2 + ,3922) ->• wJ/ijj) = 11.8 eV . (30) 

Comparing these results with the experimental predictions in Eq. ([27)1 and considering 20% 
uncertainties in the 77 rates for the (2 + , 3922) state, the results of Eq. ( 1301) are compatible 
with the assumption that the X(3915) is the resonance (2 + , 3922), considered as the Z(3930) 
in Table QUI 

On the other hand, if we assume that the X(3915) corresponds to our (0 + , 3943) reso- 
nance, the discrepancies are very large, even if we accept a factor three uncertainty in our 
result for T 77 . Our study, thus, favors the association of our (2 + , 3922) resonance to the 
X(3915) of ; 



IV. CONCLUSIONS 

We studied the two-photon and photon-vector meson decay properties of dynamically 
generated resonances from vector-vector coupled channels in a unitarity hidden gauge for- 
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malism. In the present work the focus was set on the hidden-charm mesons around 4 GeV 
analyzed in |8j. According to their masses and widths three of them are good candidates for 
the F(3940), Z(3930) and X(4160) mesons discovered by BELLE and BaBar. Further on 
the Z(3930) possibly corresponds to the recently observed X(3915). The two-photon decay 
width of the y(3940) in the hidden gauge formalism, assumed to be + , is more uncertain 
due to large cancellations, and could be compatible with the larger width of the prediction 
obtained in the pure D*D* molecule interpretation in fl7|. In the case of the X(4160), which 
we assume to be the 2 + state at 4169 MeV, the 77 decay width agrees with the results of 



171 ] . for a D*D* molecule with a chosen mass of 4140 MeV. 



The information on the T 77 decay rate of the X(3915) favors the association of this 
resonance to the (2 + , 3922) resonance that we obtain. Here, the quantum numbers + are 
clearly disfavored. 

Unfortunately there is not much data on radiative decays of the X, Y and Z mesons. 
The large variety of results obtained by us concerning the different decays and different 
resonances indicates that these measurements are very useful to shed light on the structure 
of these resonances. In particular the vector meson-photon decay modes could be addressed 
in future facilities like PANDA and BESIII, and undoubtedly these measurements would be 
very valuable to help to disentangle different aspects of the nature of these charmonium like 
states. 
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